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PART A. 

MODEL FOR REGULATION OF ERYTHROPOIESIS 



Model for Regulation of Erythropoiesis 


A mathematical model of the erythropoiesis 
regulating system was developed for use as a tool to 
investigate the relative influence of the controlling 
factors of erythropoiesis on total red blood cell 
(RBC) mass The loss of red cell mass has been a 
consistent finding during space flight and certain 
ground-based analogs of weightlessness such as bed 
rest. Computer simulation of this phenomenon re- 
quired a model that could account for oxygen 
transport, red cell production, and red cell destruc- 
tion. 

The control of erythropoiesis is amenable to a 
modeling approach, as demonstrated by a variety of 
investigators These models have ranged from con- 
ceptual. qualitative feedback schemes (refs 111-51 
and 111-52) to more de ailed mathematical descrip- 
tions (ref 111-53). some of which have led to com- 
puter simulation models (refs II1-8 and III-54). 
Quantitative models have been developed for 
different aspects of erythropoiesis regulation, in- 
cluding ferrokinetic models (refs 111-55 and III-56) 
and those describing the control of stem cell pro- 
duction (refs 111-57 to 111-59). 

A limited number of models are available that 
integrate the entire feedback circuit in sufficient 
detail to enable quantitative simulation of diverse 
hematological stresses directly affecting oxygen 
transport, tissue oxygenation, erythropoietin 
release, and red cell production. Hodgson (ref. 
111-53) has proposed such a model but apparently 
has neither implemented it for computer use nor 
performed the studies necessary for validation. An 
algorithm capable of predicting red cell mass 
changes by feedback processes was coupled to a 
much large: mode! of circulatory and fluid control 
(ref III-8), but it provided only a gross representa- 
tion of the renal/bone-marrow axis and failed to 
recognize renal oxygenation as a major control ele- 
ment The most complete simulation mode! pre- 


viously developed is that of Mylrea and Abbrecht 
(ref. lil-54); this model was validated for a single 
stress (altitude hyooxia in mice). 

A review of the factors that control 
erythropoiesis and of existing models led to the for- 
mulation of a conceptual model and eventually the 
design and development of a computer model (ref. 
111-60) The initial structure was derived from the 
best features of the mathematical models of 
Guyton et al. (ref. 1II-8), Hodgson (ref. 111-53). and 
Mylrea and Abbrecht (ref 111-54). The model has a 
wide range of applicability, the most important ap- 
plication is its simulation of space-flight events 

Model description — Those elements in the feed- 
back regulation loop that have been incorporated 
into the mode! are shown schematically in figure 
111-27 The formulation is based on the accepted 
concept that erythrocyt* production is governed by 
the balance between oxygen supply and demand in 
the body and. in particular, at a renal sensing site. 
The mechanisms and pathways of the control cir- 
cuit include oxygenation of hemoglobin and ox- 
ygenation of tissues by blood transport and diffu- 
9icnal processes An erythrocyte-stimulating hor- 
mone, erythropoietin, is produced and released at a 
renal site in amounts that vary inversely with the 
levels of tissue oxygen tension (pO:). Plasma levels 
of erythropoietin control the production of 
erythrocytes in the bone marrow. The amount of 
circulating red cell mass is based on the addition of 
new cells to existing cells with allowance for the 
destruction of older cells Blood hematocrit, an in- 
dex of red cel! concentration in whole blood, is 
determined by dividing the circulating red cell mass 
by the plasma volume. Other features of the model 
include an oxygen-hemoglobin affinity, which can 
be varied, and time delays which represent time for 
erythropoietin distribution in plasma and time for 
maturation of erythrocytes in the bone marrow. 


*This section was originally written as part of a more comprehensive chapter 
on physiological models which is as yet unpublished. Figures, tables, 
equations, numbers, and references have been left intact from the original. 
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In terms of control system theory, the controlled 
system includes the elements of redistribution of 
new cells, lung oxygenation, blood transport, and 
tissue oxygen extraction, whereas the controller 
system consists of erythropoietin release, marrow 
red cell production, and splenic destruction. The 
hematocrit (or hemoglobin level) may be con- 
sidered to be the primary feedback quantity, and 
the level of tissue oxygenation can be taken as the 
directly controlled variable. Figure 111-2“ indicates 
the input parameters which have been found most 
useful in these studies, as well as the predicted out- 
put responses which the model is capable of 
generating. Appendix B and Section V contain a 
more complete description of the model system. 

Sensitivity analysis — The effect of variation of 
the mode! parameters on the overall behavior of 
the system can be studied by sensitivity analysis 
techniques (See app. A.) Sensitivity analysis is a 
systematic method of evaluating the relative im- 
portance of different parameters. Steady-state sen- 
sitivity coefficients were estimated by varying each 
of the most important parameters, one at a time, by 
small increments around the normal control state. 
Table III-6 contains the percent change in one de- 
pendent variable (red cell production rate) caused 
by a 1 -percent decrease in the model parameters 
(shown in the left column) and determined after a 
new steady-state condition was reached Positive 
values of these coefficients indicate a response to a 
hypoxic condition, and negative values indicate 
response to tissue hyperoxia. In each case, the 
response is in a direction that agrees with the con- 
cepts of erythropoiesis regulation (ref. III-61) The 
unusually high influence of oxygen uptake on red 
blood cell production indicated may not be repre- 
sentative of the real intact system Compensatory 
action by cardiorespiratory factors usually inter- 
vene to normalize oxygen transport and thereby 
limit the correction required by the erythropoietic 
system (app. B) However, this analysis does sup- 
port the assertion that PSO shifts can dramatically 
augment tissue oxygenation (refs. III-62 and 
III-63). 

The results from a dynamic sensitivity analysis 
are illustrated in figure fII-28 In this case, the six 
most important parameters were changed by 10 
percent in a direction that would initiate a hypoxic 
response The dynamic behavior of four important 
dependent system variables during this constant 
load disturbance is shown over a period of 30 days 
It is apparent that the general ranking of parameter 
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importance found in the steady state is maintained 
for the dynamic state as well. If, however, a new 
equilibrium state was simulated for a specific hy- 
poxia condition and sensitivity coefficients were 
obtained by small perturbations from that state, 
both the magnitude and the ranking of the coeffi- 
cients might be altered. Such behavior would not be 
unexpected for a nonlinear control system. 

Figure III-28 also illustrates the wide difference 
in dynamic properties of the various elements in 
the model. This model may be viewed as consisting 
of four sequential processes, each with a charac- 
teristic time constant that is correspondingly longer 
with respect to alterations in tissue oxygenation 
(seconds to minutes), erythropoietin release 
(hours), red cell production (days), and red cell 
mass (days to months). An interesting prediction 
of these simulations is that the dynamic charac- 
teristics of the hypoxia response are similar for a 
variety of constant load disturbances, whether it be 
an increase in plasma volume or a decrease in renal 
flow rate Unfortunately, the experiments required 
to confirm this conclusion are difficult to perform. 
These simulation results should, therefore, be con- 
sidered as model-to-model comparisons which 
demonstrate that the mode! responds appropriately 
in a gross sense to both equilibrium and dynamic 
states 

Model validation— Model-to-data comparison is 
a more challenging validation process which must 
be used to establish ultimate credibility of any 


model system The erythropoiesis control model 
has been validated Tor long-term altitude hypoxia. 

Buderer and Pace (ref. 111-64) studied (he 
dynamic changes in red cell mass, hematocrit, and 
plasma volume in sea-level pigtailed monkeys dur- 
ing 6 months at 3800 meters altitude followed by 
descent to aea level for 3 months. Comparable dau 
for human subjects were not available. However, 
the human model would not be expected to be fun- 
damentally different from that for a monkey at sea 
level, especially if responses are expressed as nor- 
malized values. The major driving function for the 
model was an arterial blood oxygen tension of 6.7 
and 12.7 kN/m J (SO and 95 ramHg) for the altitude 
and sea-level phases, respectively. The model's 
response was adjusted using the bone marrow con- 
troller gain until a visual “best fit” was obtained. 
This is an example of parameter estimation. (See 
app. E.) The comparison of me del output and data 
are illustrated in figure 111-29. Iu combination with 
the arterial blood oxygen tension, the experimen- 
tally determined change in plasma volume (bot om 
curve) was used as a driver for the simulation, but 
the effect on the overall response system for this 
secondary load input was relatively small. 

In addition to predicting the measured dynamic 
changes of red cell mass and hematocrit, the model 
is capable of predicting other system variables that 
were not measured, such as plasma erythropoietin 
and red cell production-destruction rates. The 
simulation shows the general sequence of events, 
generally assumed to characterize the hypoxic 
response, including tissue hypoxia (not shown), 
elevated erythropoietin, and augmented red cell 
production levels, which increase the mass of cir- 
culating blood cells (refs. III-65 to III-67). Figu r e 
111-29 also shows the opposite scenario for the de- 
scent phase. The difference between production 
and destruction rates provides a visual indication of 
the dynamic behavior of deviation from steady 
state The slow approach to equilibrium at altitude 
is also evident from the asymptotic nature of the 
measured quantities. Upon return to sea level, the 
increased hematocrit serves as a prolonged 
stimulus for tissue hyperoxia to the extent that red 
oell production may be totally inhibited for several 
weeks. 

The erythropoietin response in humans and 
mice is known to return toward control more 
rapidly than indicated in figure III-29 (refs MI-65 
and MI-68) Adaptive mechanisms of the circulato- 
ry, ventilatory, and biochemical systems respond to 


hypoxia and thereby improve oxygen transport 
faster than does the more sluggish erythropoietic 
system (ref MI-69). Several simulations were per- 
formed to test the hypothesis that some of these 
pathways contribute to the observed erythropoietin 
response Figure MI-30 illustrates the relative in- 
fluence of three adaptive mechanisms (an increase 
in pO: due to greater ventilation efficiency, a 
decreased plasma volume, and a diminished ox- 
yhemoglobin affinity) as predicted by the model. 
These simulations, although not strictly a valida- 
tion study, not only illustrate a more realistic 
response but also indicate the manner in which 
simulation techniques can be used to test hy- 
potheses regarding physiological mechanisms. 

Mouse model — Experiments with animals are 
routinely conducted in support of space-flight 
biomedical investigations. In particular, mice have 
been examined for their potential use as experi- 
mental models for the study of the early 
erythropoietic disturbances found in astronauts 
(ref. 111-70) The desirability of using a species- 
specific model of erythropoiesis for predicting ex- 
perimental results with animal subjects prompted 
the formulation of a “mouse model." This mode! 
was based on the human model but incorporated 
parameter values that describe the unique charac- 
teristics of the oxygen transport and erythropoietic 
system for the mouse Thus, red cell lifespan, ar- 
terial oxygen tension, oxygen-hemoglobin affinity, 
and many other parameters were found to differ 
significantly from the human (ref. MI-71 and app. 
B). 

The mouse model was validated in a preliminary 
study by examining its responses to hypoxia (ref. 
III-70) and to erythropoietin infusions (ref. IM-72). 
Ar example of the dynamic behavior of the mouse 
model is provided in figure III-31. In this simula- 
tion, a pulse decrease in arterial oxygen tension 
results in a conesponding disturbance in tissue ox- 
ygen levels, followed by a short burst of 
erythropoietin and t resulting wave of 
erythropoiesis, the effects of which may last 2 
weeks or longer The transient nature of red cell 
production shown here has been previously ob- 
served after injecting erythropoietin into mice (ref 
III-73). A similar simulation performed with the 
human model results in more prolonged effects, as 
indicated in figure 111-31 These ifTerences result 
primarily from the shorter erythropoietin disap- 
pearance rates, red cell lifespan, and erythropoietin 
half-life in the mouse when compared to the 
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human. Simulations such as this suggest that some 
information regarding erythropoiesis may be ob- 
tained in a shorter time interval using the mouse as 
an experimental model of the human 

Other simulations — Although the major purpose 
of the erythropoietic control model was to investi- 
gate mechanisms related tc the loss of red cells dur- 
ing space flight, other applications have been 
investigated to validate the model, to test hy- 
potheses, and to introduce refinements Table III-7 
is a summary of the major simulations ac- 
complished with the human and mouse models 
Human bed-rest and mouse dehydration studies 
have been found ;o be useful experimental analogs 
of space-flight efTects. A more recent area of model 
development is in the simulation of hematological 
disease states which have known physiological 
etiologies (such as anemias, polycythemias, and ab- 
normalities in hemoglobin functions) Toward this 
end, a steady-state model based on the dynamic 
model described previously has been formulated A 
program is now in progress to develop a teaching 


model of pathophysiology based on the human red 
cell model. 

Limitations — The erythropoiesis control system 
model is the least complex of all the models de- 
scribed in this section. Various features have been 
added during its development to increase its simu- 
lated realism These have included elements repre- 
senting erythropoietin release, cell maturation time 
delays, and a variable oxygen-hemoglobin affinity 
Other limitations have been identified which either 
are acceptable within the scope of this program or 
form the basis for suggested improvements These 
limitations are summarized here. 

1. Red cell production is modulated by factors 
other than those considered in the model These in- 
clude iron levels, hormones other than 
erythropoietin, and neural stimulation, as well as 
inhibitors and activators of erythropoietin 
Although it may be desirable to include these 
effects in future model applications, quantitative 
information regarding them is presently lacking or 
the need for their inclusion is not yet warranted 
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2. The model has a limited representation or 
bone marrow red cell production For example, it 
does not include the full maturation cycle from 
multipotential cells to denudeated red cells 
Furthermore, it has no means of predicting 
reticulocyte index. 

3. The model does not include extrarenal 
erythropoietin production or extramedullary 
erythropoiesis These factors are of much greater 
importance in the mouse than in the human 
However, sufficient quantitative information is 
lacking for inclusion in a mathematical model 

4 Blood viscosit> and blood volume changes 
are known to influence blood flow and. thereby. 


tissue oxygenation These effects are not explicitly 
represented in the model (refs. 111-81 and 111-82). 

5 Overall oxygen transport is under the control 
of homeostatic mechanisms in addition to the 
erythropoiesis system Thus, some of the fixed 
parameters of the model such as blood flow, capil- 
lary difTusi vity. arterial oxygen tension, p'asma 
volume, and oxygen-hemoglobin affinity can be 
considered variable elements of circulatory, ven- 
tilatory, biochemical, and fluid regulatory feedback 
mechanisms that are beyond the scope of the pres- 
ent models design objectives. 

Larger models that incorporate many of these 
features are available, and these have been shown 
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Table III - 7.- -Computer Simulation Studies 
Performed With Erythropolesls Model 


Study 

Duration. 

days 

Scurcr 

Human model 


Altitude hypoxia and descent 

270 

Ref 111-64 

Altitude nyponia 

6 

Ref 111-67 

Red cell infusion 

120 

Ref HI-74 

Bedrest 

35 

Rers 111-75. 
111-76 


14.28 

Ref HI-77 

Skylab 

59.84 

Ref 111-78 

Clinical conditions 


” 

Anerva (primary, sec- 
ondary, hemolytic) 
Polycythemia (primary. 

secondary, hypoxic) 
Abnormality (hemoglobin) 

Mouse model 


Multiple sequential 

erythropoietin 

7 

Ref 111-73 

infusion 

Altitude hypoxia 

28 

Ref UI-bJ 

Dehydration and recovery 

7 

Ref 111-79 

Red cell infusion 

3 

Ref III-80 
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1979) 
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PART B. 

DETAILED DESCRIPTION OF THE MODEL FOR REGULATION 
OF ER YTHROPO I ES I S 


* Detailed Description of the Model for Regulation of Erythropoiesis 


The model developed for the space-flight 
hematological program was designed to examine 
the relative influence of the controlling factors of 
erythropoiesis on total red blood cell (RBC) mass. 
The formulation was based on the generally ac- 
cepted concept that the overall balance between ox- 
ygen supply and demand regulates the release of a 
hormone, erythropoietin, from renal tissues sensi- 
tive to oxygen tension levels and that this sub- 
stance in turn controls bone marrow red cell pro- 
duction (refs. B-l to B-3). 

The amount of oxygen delivered to the tissue is 
accounted for in the model by the combined in- 
fluence of several factors: hemoglobin concentra- 
tion, lung oxygenation of hemoglobin, blood flow, 
and oxygen-hemoglobin affinity (ref. B-4). From 
this amount of oxygen, a certain fraction is ex- 
tracted by the tissue, depending on the oxygen de- 
mand parameter. Oxygen enters the cellular spaces 
by diffusion along an oxygen tension gradient be- 
tween the venous capillaries and the cell (ref. B-5). 
Decreasing the oxygen supply in relation to the de- 
mand reduces tissue ox: gen tension which is, in 
effect, monitored by a local oxygen detector (ref. 
B-6) and results in increased rates of erythropoietin 
release (ref. B-7). Erythropoietin is released into 
the general circulation, and its final plasma con- 
centration is determined by its rate of release, its 
volume of distribution, and the rate at which it is 
metabolized (ref. B-8), the latter being represented 
in the model by the hormone piasma half-life. The 
target organ for erythropoietin is the bone marrow. 
The production rate and release of red cells in the 
model are determined by the plasma erythropoietin 
concentration (ref. B-9). A time delay exists be- 
tween marrow stimulation and red cell release (ref. 
L-10). Destruction of cells is represented in the 
model by a lifespan parameter. Hemoglobin con- 
centration in blood is based on the addition of new 
cells to existing cells and plasma. 

It is known that red cell production is modulated 
by factors other than those considered in the 
model. These include iron levels (ref. B-l I), hor- 
mones other than erythropoietin (ref B-l 2), neural 
stimulai (ref B-l 3), as well as inhibitors and ac- 
tivators of erythropoietin (ref B-14). Although it 
may be desirable to : -Hude the. effects in future 


model applications, quantitative information 
regarding them is presently lacking or the need for 
their inclusion is not yet warranted. 

Overall oxygen transport is under the control of 
homeostatic mechanisms in addition to the 
erythropoiesis system. Thus, some of the fixed 
parameters of the model such as blood flow, capil- 
lary djflusivity, arterial oxygen tension, plasma 
volume, and oxygen-hemoglobin affinity can be 
considered variable elements of circulatory, ven- 
tilatory, biochemical, and fluid regulatory feedback 
mechanisms that are beyond the scope of the cur- 
rent model’s design objectives Larger models that 
incorporate many of these features are available 
and these have been shown to be compatible with 
an erythropoiesis subsystem model (ref. B-l 5). 
Although these mechanisms are not included ex- 
plicitly in the present erythropoiesis model, their 
influence can b: tested, in most cases, by manual 
alteration of existing model parameters. 


MATHEMATICAL DESCRIPTION OF MODEL 
Oxygenation of Blood 

The oxygen concentration of arterial blood 
(neglecting oxygen dissolved in plasma) after 
passage through the lungs can be expressed as the 
product of the carrying capacity of a gram of 
hemoglobin ( CHbO ; i.e., 4 moles 0: per mole Hb), 
times the hemoglobin concentration ( Mb ), times 
the fractional degree of oxygen-hemoglobin si tura- 
tion (S fl 0). It is convenient to use hematocrit rather 
than hemoglobin concentration as an index of red 
cell concentration in blood. These quantities are 
related by the mean corpuscular hemoglobin con- 
centration, such that MCHC - Hb/Hct. Therefore, 
for arterial blood 1 

C'O * S # 0 X Hb X CHbO 

* S'O X Hct X MCflC X CHbO (Bl) 


'Set table D-l Tor symbol definitions, units, and normal 
values 


♦This section was originally written as part of a more comprehensive chapter on 
physiological models which is as yet unpublished. Figures, tables, equations, 
numbers, and references have been left intact from the original. 


B-l 




At full hemoglobin saturation, 5,0 m 1 .0 and equa- 
tion (Bl) then represents the maximum 
oxygen-carrying capacity of blood at a given 
hematocrit. In most instances, the parameters 
A/C//Cand CHbO will be invariant, so that the ar- 
terial oxygen concentration is influenced only by 
the hematocrit and hemoglobin saturation In prac- 
tice, the arterial oxygen tension (P a O) is assigned a 
value, 1 and a corresponding value of 5,0 is deter- 
mined from the oxygen-hemoglobin equilibrium 
curve (0£Q, an operation tha: can be expressed in 
functional form as 

S'O - OEC(P g O) (B2) 


Oxygon Delivery at Tieeuea 

Oxygen transport at the tissue level is schemat- 
ically represented in figure B-l. In the present 
model, the tissue of concern is taken to be the 
erythropoietin-producing cells known to be located 
primarily within the kidneys It has been assumed 
thM the oxygen sensor as well as sites of 
erythropoietin production are responsive to the 
mean oxygen tension of the kidneys This assump- 
tion permits tissue oxygen tension to be derived 
from an oxygen balance using blood flows, ar- 
teriovenous oxygen concentrations, oxygen con- 
sumption, and transcapillary diffusion resistances 
common to the kidneys, which are considered to be 
a homogenous tissue. 

'See ubte B-l for symbol definitions, units, and normal 
values 



V m 


FIGURE B- 1 . — if rtMcn t Mredt-stair aedel ef renal 
Itaeve n?|mllw 


The elements of tissue oxygenation were deter- 
mined using a two-phase model consisting of capil- 
lary blood and tissue fluid (ref. B-16) It is assumed 
that the blood compartment is well mixed with an 
oxygen partial pressure ( P ¥ 0) equal to that in the 
venous outflow Oxygen diffuses from blood to 
tissue along a gradient of oxygen partial pressute 
(PjO - P t O ), where P,0 is the oxygen tension of 
the homogenous tissue. The amount of oxygen 
unloaded from the blood is given by the difference 
in oxygen concentration between arterial and 
venous blood (C,0 — CjO ). At equilibrium, the 
rate of oxygen transfer to the tissues is identical to 
the tissue oxygen consumption ( VJ and to the un- 
loading of blood oxygen as determined by the ar- 
teriovenous oxygen concentration In other words, 

Q x (C g O - C v O) = K d X ( P f O - P,0) 

Rate of oxygen Rate of oxygei 

unloading diffusion 

* Kn ™ 

where Qis the rite of regional blood flow and K d is 
the diffusive transfer coefficient between blood 
and tissue. 

The venous partial pressure may be obtained by 
solving equation (B3) for C v O 

C,0 - C,0 - (VJQ) (B4) 

then computing venous hemoglobin saturation 

from a formulation analogous to equation (Bl) 

S f O « C,0/(Hct X MCHC X CHbO ) (B5) 

and determining PJO from the OEC, expressed in 
functional form as 

PJO = OEC(S r O) (B6) 

Tissue oxygen tension can now be derived from 
equation (B3) as 

ffi - f,0 - V m jK , <B7) 

The assumption of steady state in this formula- 
tion implies rapid equilibration of oxygen in the 
fluid phase and is not meant to suggest a constant 
tissue oxygen tension during the simulation of an 
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erythropoietic disturbance Since the erythropoietic 
process itself has been modeled dynamically, the 
hematocrit will be time-varying until red cell pro- 
duction achieves its own steady state In this 
model, the hematocrit is a major influence on 
alterations in tissue oxygen tension In addition, 
other quantities in this algorithm (PjO, V m , 

£50, Q) are nonregulatory parameters and can also 
influence tissue oxygen tension if they are 
manually altered. 

Erythropoietin Production and Distribution 

Erythropoietin release is assumed to be 
governed by the tissue oxygen tension. A limited 
number of studies (refs. B-5 and B-17 to B-19) sug- 
gest a relationship of the form 1 

r p • (B8> 

where £Jis the rate of erythropoietin production, 
P,0 is the tissue oxygen partial pressure, Gi is the 
gain or slope ofthe function plotted linearly as the 
natural log of E p versus P,0 , and E 0 is the ^inter- 
cept at Pp - 0 Setting E 0 - exp(Gi) ensures that 
equation (B8) wiljjlways pass through the normal 
operating point (E p - 1.0, P,0 - 10) irrespective 
of values of G In certain cases, it may be desirable 
to postulate a shift in the norma! operating point as 
well as gain, and in those cases, E 0 and G may be 
adjusted separately. Figure B-2 illustrates the semi- 
logarithmic relationship between tissue oxygen ten- 
sion and erythropoietin production and shows the 
influence of G. 

The concentration of erythropoietin in the 
plasma (£) is a function of the rate of production 
(£,), the rate of clearance or destruction (EJ, and 
the volume of distribution (F r ). If it is assumed 
that the rate of disappearance is proportional to the 
plasma concentration (i.e., E d — K t x V f x £. 
where K t - clearance constant - log, 2/plasma 
half-life), then the following first-order differential 
equation can be written for the rate of change of 
erythropoietin concentration (refs B-8 and B-17): 

(B9) 

dt V f * 

'The, bar over ihe lymboti represents normalized values 
Thus. X - JT/JTto). where X{o) is ihe control or pi«!ircaimen: 
condition and X is » value during the paciconuol or treatment 
phase The use of normalized values simplified comparison of 
model response wuh data from various laboratories 



FIGURE B-2. — Renal controller (unction for erythropoietin 
production rate showing effect of gain G. 


The steady-state concentration of £ (at dEldt • 0) 
is, therefore 

E{p) = E p {o)/K'V t (BIO) 

Equation (B9) can be normalized by letting 
£ ™ £/£(o) and E p “ E p /EJo) and substituting 
equation (Bl) into equation (B9): 

# ■ - c ) 

(b,i > 

a form that has the advantage of being independent 
of distribution volume and, further, independent of 

erythropoietin half-life ( TE xn ) at steady state, 

since £(o) - E p (o). The value of TE xn — 12 hours 
chosen for humans, but this does not appear to 
be well established for normal subjects (ref B-20) 

Rod Blood Coll Production 

An equation expressing the semilogarithmic 
dose-response relationship between red cell pro- 
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duction (RCP) and erythropoietin concentration 
(£) can be given as 

RCP = C 2 Iob, E ♦ />, (B12) 

As has been the convention, RCP and £h«ve been 
normalized with respect to their steady-state con- 
trol values, Gi is the gain or slop e of the response 
curve, and P\ is the value of RCP when £“1.0 and 
is normally taken as unity. 

The accuracy of equation (B12) may diminish 
considerably at very low and very high values of 
erythropoietin levels. For example, at decreasing 
values cf £ approaching zero, the production rate 
tends toward minus infinity rather than zero, 
whereas at high £ values, the relationship does not 
exhibit a maximal production that is known to ex- 
ist. Therefore, two additional equations were for- 
mulated and piecewise fitted to equation (B12) to 
account for suppressed and maximal erythro- 
poiesis (See fig. B-3.) Their precise description r 
somewhat speculative, although this bone marrow 
function curve closely corresponds to the sigmoid 
shape of most biological dose-response relation- 
ships (ref. B-21) and to those observed by others 
for the bone marrow (refs B-22 to B-24) 

In the present study, it was assumed that there is 
no basal production of red cells unless erythro- 
poietin is present (i.e., P 0 “ 0). Maximal red cell 
production is assumed to be six times normal (ref. 
B-17) (P m “ 6) and the upper limit of ac curac y of 
equation (B12) was arbitrarily chosen at RCP “ 5. 

Iron uptake studies of red cells indicate a bone 
marrow transit time for red cell production of 3.S to 
4.5 days (ref. B-10) This effect was included in the 
model using a simpie first-order time delay 1 with a 
time constant Tg M . The inclusion of this transit 
delay clearly improved ' he realism of the dynamic 
response of the model, especially following sudden 
changes in tissue oxygen tension. 

Rpd Cell Destruction and Rod Coll KAsts 

The lifespan of the red cell dicutes the destruc- 
tion rate of cells A model of red cell destruction 


'If Y a iht steady -m tc value predicted from a dose 
response relationship. th»n the response delayed with a time 
constant T is ohtair^-d o> the solution of the difTerennal equa- 
tion T dwdi 4i« } In finite difference form, suitable for itera- 
tive computer solution, this becomes jn m jn + ( Y - 
where H - inicfration step size. * value of i at the ah itera 
non. and - value of » at the i 4 Ith iteration 


was assumed in which cells are destroyed randomly 
and no cells survive past a given lifespan (ref. 
B-25) In that case, the rate of destruction (RCD) is 
simply proportional to the amount of cells present 
at any time If that amount is given as RCM , then 
the rate of red blood cell destruction will be 

RCD * K f ‘RCM 

= log, 2 -RCM/TRC^ (B13) 

where TRC l/2 “ red cell half-life and K, “ red cell 
clearance consunt — log, 2 JTRC m . At normal 
values used in the model, the rate of red cell 
destruction is 1 1 percent of the toul amount pres- 
ent or 22 ml of packed cells/day. 

The instantaneous change in toul circulating red 
cell mass is the net difference between red cell pro- 
duction and red cell destruction rates: 

= RCP - RCD (B14) 

at 

The quantity RCP is obtained fro m equ ation (B12) 
using the transformation RCP “ RCP x RCP{o ), 
where RCP[o) “ RCD(o ) - K,-RCM(o). 

Finally, the current value of RCM , obuined by 
integration of equation (B14), is combined with the 
plasma volume to obtain the whole-body 
hematocrit of the circulation: 

*• ■ wrm i (B15) 

and the feedback loop is closed. 


Th# Renal/Bonc-Marrow Axis at a 
Proportional Controflar 

A more visual understanding of the relation- 
ships for renal and bone marrow function can be 
obuined by determining the overall dose-response 
curve of the combined kidney/bon e-marrow axis; 
that is, the relationship between tissue oxygen ten- 
sion and red cell production. The results are shown 
here only for the midrange of erythropoietin pro- 
duction and were obuined by mathematically com- 
bining the renal function (eq (B9)) with the 
steady-state plasma distribution function of 
erythropoietin (eq (Bll), £ - £^) and the bone 
marrow function (eq (B12)). In this way. 
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FIGURE B-3. — B*nf Barra* c*o trailer faoction far red cell production rale showing equations used to p i e cew is e fit the dose- 
response curve and demo ut rating effect of gain ft. 


erythropoietin is eliminated as an explicit variable; 
i.e., 

RCP = G(l - Pfi) (fot Pfi < l) (B16) 

where G ™ Gi x Gi. This is a simple inverse linear 
relationship that has previously been used in the 
models of Guyton et al. (ref. B-l 5) and Hodgson 
(ref. B-5). In this form, the gain or the renal/bone- 
marrow (open loop) system is seer to be the prod- 
uct of the two gain factors, G and G> Similar 
(although nonliuear) gain-product relationships 
were obtained for the extreme of hyperoxia and 
maximal hypoxia. The results of this analysis have 
been plotted in figure B-4 for a range of values of G 
This function represents the entire controller cir- 
cuit of the erythropoiesis model. In general, hypox- 
ic tissue is predicted to have a greater effect on red 
cell production than hyperoxic tissue, as indicated 
by the relative slopes. However, at higher gain fac- 
tors, the suppression of erythropoiesis can be con- 
siderable for relatively small increases in tissue ox- 
ygenation. 

Equation (BI6) is a typical controller function 
for a linear proportional control system in which 
the actuating error signal is the deviation of P,0 
from its control value (unity) and G is the constant 
of proportionality or controller pm The control 
value of Pfi is taken as an arbitrary reference 
standard and is not meant to imply that there is an 
internal set-point It becomes clear from this con- 


trol function that the primary controlled variable is 
tissue oxygen tension rather than red cell produc- 
tion, red cell mass, or hematocrit, which can be 
considered controller variables. The control system 
always adjusts the controller variables in a direc- 
tion that returns the controlled variabie toward 
normal after an initial load disturbance (See fig. 
A-l in app. A.) 

The controller gain is a major determinant of the 
final feedback compensation following a dis- 
turbance and, within limits, the speed at which a 
disturbance is corrected Although increasing the 
pin increases the feedback effect, it may be at the 
sacrifice of an oscillatory approach to the steady 
state Systems of higher order than this one will 
tend to become unstable as controller gain in- 
creases Oscillations in the erythropoiesis system 
are known to occur in special cases, nd analysis of 
these systems in terms of control system theory 
has proved to be rewarding (refs B-26 and B-27). 


Modal Operation 

The system of equations (Bl), (B2), (B4) to 
(B8), and (BU) to (BIS) is solved by computer 
simulation using an iterative procedure. Equations 
(B’.l) and (BI4) are integrated numerically using a 
simple Euler algorithm (ref B-28) and initial condi- 
tions E(o) and RCM(o), respecrively The program 
is currently implemented ofi a Univac 1110 and a 
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FIGl'RE B-4 —Combined renal /bone-marro« controller function showing effect of tissue oxygen tension on erythropoiesis for 
various values of overall gain G. 


PDP 11/40 using Fortran language Remote ter- 
minals including graphic display of the simulation 
responses (and experimental data simultaneously, 
when desired) have greatly facilitated user interac- 
tion and have enhanced the convenience of model 
validation 

Perturbing the model away from its normal 
steady state (initial conditions) is accomplished by 
altering one of the model parameters. 4 These quan- 
tities are normally constant in value but can be 
altered during a run to simulate either independent 
stress stimuli or long-term regulatory adjustments. 
Thus, hypoxia may be simulated by decreasing ar- 
terial oxygen partial pressure (pO?),and hemolysis 


The term “piumeier” is used 10 denote iime-inv*ri»nt 
quinnnes the ire not altered b> the dynamic properties of the 
model as opposed to "dependent variables" thai are lime- vary- 
ing and form the connecting links of the feedback circuit lie . 
P,0 L RCP Ha etc I 


may be simulated by decreasing the red cell 
lifespan In some cases, alteration of more than one 
parameter may be desired as in the simulation of 
altitude hypoxia which involves a primary change 
in arterial pO: and a compensatory shift in oxygen- 
hemoglobin affinity (P50). Parameter values may 
be time-invariant with values fixed before a simula- 
tion run or they may be entered as a function that 
varies with time as the run progresses Experimen- 
tal data may be used to drive the model in this 
fashion Thus, it is possible to simulate a wide 
variety of stresses and to test a large number of hy- 
potheses regarding regulating mechanisms. 


Parameter Estimation 

Table B-I lists the values of the system 
parameters, chosen as representative of the human 
system Also shown are the steady-state control 
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Table B-l. — Parameters and Initial Conditions 


Symbol 

Definition 

Value 

inns 


Parameters 


CHbO 

Carrying capacity of hemoglobin 

1.34 

ml 0,/g Hb 

f. 

Intercept of erythropoietin function 

MW 

x normal 

c, 

Cam of erythropoietin function 

3 

Nondimensional 

G, 

Gain of red cell production function 

2 

Nondimensional 

*4 

Capillary diffusivity 

.567 

mlO,min~ l mmHg* 1 



425 

tnlOimin^'kPa -1 

MCHC 

Mean corpuscular hemoglobin concentration 

.375 

g Hb/ml RBC 

P.0 

Oxygen tension in arterial blood 

950 

mmHg 



12 67 

It Pa 

PV 

Plasma volume 

3.0 

Liters 

P m 

Maximum production rate of red cells 

6 

x normal 

P. 

Minimum production rat. of red cells 

0 

x normal 

P, 

Normal production rate of red cells 

1 

x normal 

Pi 0 

Oxygen tension of hemoglobin at SO percent saturation 

26 7 

mmHg 



356 

kPa 

Q 

Rena: blood flo» 

12 

Liiers/min 

TBM 

Bone marro* transit time 

4 

Days 

TE,r. 

Plasma half-life of erythropoietin 

12 

Hours 

TRC xn 

Red ce" half-life 

63 

Days 


Oxygen upuke of kidneys 

M 

ml 0,/min 

Initial conditions 

C.O 

Oxygen content of arterial blood 

196 

ml 0,/l 00 ml blood 


Oxygen content of venous blood 

|7 9 

ml O.MOO ml blood 

£ 

Erythropoietin concentration in plasma 

1.0 

x normal 


Erythropoietin production rate 

1.0 

x normal 

Ha 

Hematocrit 

40 

ml of packed RBC'100 ml 




blood 

PjO 

Oxygen tension ir. renal tissue fluid 

2C 

mmHg 



2 67 

kPa 

P.0 

Oxygen tension in venous blood 

55 3 

mmHg 



7.37 

kPa 

RCM 

Red cell mass 

20 

Liters 

RCP 

Production rate of new red blood cells 


ml of packed cells/day 

S.0 

Saturation of arterial hemoglobin with oxygen 

97 6 

Percent 

s.o 

Saturation ei venous hemoglobin with oxygen 

893 

Percent 


values (i.e., initial conditions) of the major de- 
pendent (output) variables The precise values of 
most of the quantities shown in table 6-1 are not 
critical, however, to the general behavior of the 
model’s response when expressed in percent devia- 
tion from control An exception to this occurs for 
the highly nonlinear functions such as the oxygen 
equilibrium curves and the renal and bone marrow 
function curves The location of the operating point 
on these curves can change the nature of the output 
response Also, the shapes of these curves, as 
defined by the gains and threshold parameters are 
critical to '.he magnitude of the response These will 


be discussed next. 

An equation for the oxygen equilibrium curve 
was obtained from Aberman et al (ref. B-29). Their 
study included a computer algorithm that converts 
oxygen tension to saturation and, with iteration, 
oxygen saturation to tension The agreement be- 
tween measured and computed values is claimed to 
be better than 0 2 percent In addition, the 
algorithm includes the capability of altering P50 
values. (See fig B-5 ) 

Analysis of Adamson's data (ref B-l 8), from 
which the erythropoietin release function was ob- 
tained, leads to an estimated value of G " 2 8 In 
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pOj, mmHg 

I 1 1 1 1 1 

0 2.67 5.33 8.0 10.67 13.33 

i0 2 . kN/.* ? 

FIGURE B-5 — Oxygen-hemoglobin equilibrium runes drawn 
direct!) from computer model output showing effect of chan|- 
in| the KB parameter. The A50 value is the oxygen tension at 
which SO percent of the hemoglobin is saturated and is. 
therefore, an index of oxyhemoglobin affinity. The position of 
the OEC may vary in accordance with known effects of hy- 
drogen ion. COi. temperature, or Z.S-DPG (diphospboglycer- 
atei. The effects of DPG are of interest because of its possible 
contribution to erythrokinetie changes during bed rest or 
weightlessness. 


simulation studies of bed rest and hypoxia (refs 
B-30 and B-31), values of Gr. were determined by 
parameter estimation to range from 2 to 5. No data 
exist to confirm these estimates directly. It is per- 
missible to use the gain factors as adjustable 
parameters, within reasonable limits, in order to 
test hypotheses and possibly provide improved 
agreement between model and experiment data In 
addition to the gains or sensitivities, the other 
parameters of these functions P 0 . P\) that can 
be thought of as threshold indices are also not well 
known and are available for parameter estimation 
during a simulation Changing these latter quan- 
tities implies a shift of the normal operating points 
((£,- I) for fig b-2ARCP- 1,£- l) 

for fig B-3) The value of K d was obtained from 
equation (B7) by dividing the oxygen uptake by 
P % 0 - P,0 The value <f P,0 was arbitrarily 
assumed to be 2 67 kN/m- (20 mmHg), and P t O 
was computed from equations (B1 ). (B2), and (B4) 
to (B6) 


Steady-State Errors 

It is an inherent property of most biological 
homeostatic mechanisms in general and propor- 
tional coniro' systems in particular that thete will 
be at least some residual steady-state deviation 
from the normal operating point when the system 
is disturbed by a constant load This “steady-state 
error” will vary in size depending on the gain of the 
system and the magnitude of the load disturbance 
The sensitivity coefficients of table III-6, for exam- 
ple, indicate the steady-state eriors of red cell pro- 
duction resulting from small changes in parameter 
disturbances 

A basis for underfunding steady-state errors in 
the current mode! is as follows At any equilibrium 
state, the term dRCMIdt in equation (B14) must be 
identically zero Therefore, equations (B13) and 
(B14) imply that, at steady state, daily red cell pro- 
duction and destruction rates are equal Further- 
more, if red cell lifespan is constant, a steady-state 
alteration in red cell mass implies a proportionate 
change in red cell production. Most load dis- 
turbances are accompanied by changes in circulat- 
ing red cel! mass (fig B-l) This implies that red 
cell production and one or more of the factors that 
affect red cel! production (including tissue oxygen 
tension, erythropoietin, hematocrit, controller gam, 
and set-point) must have incurred a steady-state er- 
ror, however small. It can be appreciated from 
figure B-l that steady-state errors in tissue oxygen 
tension are reduced at the expense of larger errors 
in red cell mass and hematocrit. 

Steady-state analysis may also provide insight 
into the independent determinants of the con- 
trolled variable (ref B-32). Figure B-6(a) shows the 
effect of tissue oxygen tension on red cell produc- 
tion (the controller function of fig. B-4) and the bal- 
ance between red cell production and destruction 
Feedback control ensures that the various elements 
of the model are adjusted until the rate of red cell 
mass change, dRCMIdt, becomes zero and that 
steady state is defined by this condition It is math- 
ematically permissible to set dRCMIdt to equal zero 
and work backward from this point to visualize 
those factors that determine tissue oxygen tension 
(fig B-6(b)) There are only two such farors in our 
model red cell destruction (which determines pro- 
duction rate) and the relationship between tissue 
oxygen tension and red cell production Only these 
factors and those quantities that affect these factors 
will ultimately determine tissue oxygen tension 
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FIGURE B-4 Stearfy-atatr analysis of lna|-<erm determi- 

nants of tisane esyicn tension, (a) Seftnent of feedback circuit 
that results in the term itCMIii ketamine zero at stead> state 
when destruction and production rales arc equal, (b) Reversal 
of diagram la figure R-d(a) start log from the sere point and 
showing mathematical!) equivalent Influences on tissue ox- 
ygen tension. 


Experimental evidence is severel> lacking on the 
determinants of the shape of the controller func- 
tion curve, although it appears that it may be under 
the influence of neural and biochemical factors. 
fSee subsection entitled “Discussion. "j The factors 
that are responsible for normal red cell destruction 
are also not apparent (ref B-10) Except for overt 
pathology, rates of destruction are generally con- 
sidered to be a constant fraction of total circulating 
red cell mass 

Aside from these experimental shortcomings, it 
is important to note that the factors usually agreed 
to play a role in acute changes of tissue oxygenation 
(i.e., blood flow, capillary diffusivity, oxygen up- 
take, /*50, hematocrit) are not represented between 
the zero point and the tissue oxygen point in figure 
B-6 Although they may be considered dependent 
variables in the system, none of them, from a 
mathematical point of view, are independent deter 
minants of the final level at which the tissue ox- 
ygen tension will stabilize in the steady state 
Although dRCM/di will always return to its initial 
zero value, the tissue oxygen tension may exhibit a 
steady-state error due to the inherent properties of 
the control system in the face of a constant dis- 
turbance However, based on the preceding 
analysis, the authors believe that even tissue ox- 
ygen tension will return to its initial value (zero 
steady-state error) when both of the following con- 
ditions are satisfied ( I ) a constant daily rate of red 
cell destruction and (2) a constant controller rela- 


tionship between ussue oxygen tension and red cell 
production 

As an example of the last point, a simulation of 
hypoxia was performed (fig B-7) in which arterial 
oxygen tension was set at some low value for the 
entire run and red cell destruction was clamped at 
its control value of 22 ml of packed cells/day 
Tissue oxygen tension decreased and then began to 
return toward normal as red cell production and 
hematocrit rose, similar to the simulations with the 
intact system shown in figure B-7(a). However, 
since the destruction rate was not permitted to in- 
crease, a greater net rate of red ceils entered the cir- 
culation than would have occurred had destruction 
rate rose in accord with the mass action law of 
equation (B1J). As a result, tissue oxygen tension 
continued to rise and red cell production rates 
declined When the system reached its new steady 
state, red cell production and tissue oxygen tension 
returned exactly to the prehypoxic control condi- 
tions in accord with the concepts discussed in the 
previous paragraph This was despite an arterial ox- 
ygen tension that was still significantly depressed 
and at the expense of hematocrit and circulating 
red cell mass that were considerably above normal 

It is perhaps easy to visualize that a primary 
change in destruction rate, as in hemolytic anemia, 
leads to secondary changes in tissue oxygen ten- 
sion It is more difficult to conceive of destruction 
rate being a determining factor of tissue oxygena- 
tion in a stress like hypoxia in which it appears, at 
first, that the decreased oxygen loading of arterial 
hemoglobin is the primary stimulus for hypoxia 
However, it is important to distinguish between the 
initial stimulus of the acute phase, which blood pO: 
is controlling, and the ultimate stimulus of the 
steady-state condition, which destruction rate (and 
controller function) is controlling These conclu- 
sions from a theoretical model may warrant furthci 
experimental examination. The simulation of 
hemolytic anemia presented in Section VIII illus- 
trates the capability of the control system to 
minimize steady-state errors 


DISCUSSION 
Renal pO; Senior 

A major assumption in the controlled system is 
the description of the renal oxygen detector The 
evidence strongly indicates that the balance be- 
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FIGURE B-7. — Effect of destruction rule on lont-term control of tlxsue oxyien tension (a) Normal simulation response to hypoxic 
stress, (b) Hypoxic simulation with destruction rate clamped at control >aloe and showing regulation of (issue pO: bock to normal 
despite reduced f JO. 


tween oxygen supply and demand at intrarenal 
sites is the primary stimulus for erythropoietin 
release (refs B-7 and B-33). Furthermore, these 
detector sites must monitor venous or tissue pO: 
rather than arterial pO: since anemia or increased 
oxygen-hemoglobin affinity lead to increased 
erythropoietin production without significantly 
altering arterial blood oxygen tension (refs B-5, 
B-34, and B-35) The receptors sensitive to tissue 
pOj may be those cells that excrete erythropoietin 
or its precursor (ref B-13). The present model is in 
accord with these concepts. 

It has been suggested that the kidney has unique 
characterises that enable it to function as a sensi- 
tive oxygen chemoreceptor and. in particular, be 
responsive to changing hemoglobin levels The 
peculiar renal microcirculation and the uniquely 
low arteriovenous oxygen difference provide a gra- 
dient of tissue oxygen tension that amplifies 
changes in blood oxygen delivery (refs B-14 and 
B-35) In addition, the autoregulatory features of 
the kidney ensure that blood flow and oxygen up- 
take are effectively stable over a wide range of ox- 
ygen tensions and blood pressures (refs B-36 and 
B-37) Moreover, if blood flows should be altered. 


the kidney , in contrast to other organs, will exhibit 
proportionate changes in oxygen uptake (ref B-38) 
This means that the ratio I'JQ (the only term in 
which blood flow appears in the model (eq (B4)) 
may be relatively constant, and renal blood flow 
would not be expected to markedly influence tissue 
oxygenation 

The preceding discussion suggests that the 
powerful influence which oxygen uptake, per se, 
was found to have in the model (fig. B-l and table 
B-I) may be mitigated in the real system by concur- 
rent changes in blood How. Similarly, it is possible 
that an elevation of V ^ by whatever cause, pro- 
motes tissue hypoxia and results in local regulatory 
increases of the effective capillary diffusivity. 
This is now known to be true for skeletal muscle 
(ref B-39) but has not been confirmed for the 
kidney Such regulation would, however, further 
dampen the effect of oxygen uptake because K d ap- 
pears only in the ratio VjK d in equation (B7) It 
may be desirable to add these local regulatory 
effects — between oxygen uptake, blood flow, and 
capillary diffusivity — to the model A current alter- 
nate approach is to assume they are constant and to 
examine their influence, if the data so suggest, in 
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improving the accuracy of simulation Under these 
circumstances, and in accord with equations (Bl) 
to (B7), tissue pOr would be a function of the 
hemoglobin concentration of the blood, the arterial 
oxygen saturation of hemoglobin, and the shape 
and position of the OEC (ref. B-6) 

It should be emphasized that the site of the in- 
traienal detector has not yet been confirmed and 
that the quantitative aspects of its oxygen supply- 
demand balance (including direct measurement of 
P t O) remain unknown. Justification of this seg- 
ment of the model is based on indirect evidence, 
gross characteristics of the kidney as a whole, and 
determinants of tissue oxygenation derived from 
other tissues 

The use of a steady-state formulation for tissue 
oxygenation (eqs. (Bl) to (B7)) in a dynamic 
model is justified because, in the well-perfused 
kidney, equilibrium of oxygen tension due to pure 
convection and diffusion may be achieved in the 
order of seconds to minutes following a load dis- 
turbance This can bo compared to the much slower 
changes of the erythropoietin distribution or bone 
marrow red cell production process Estimates of 
true equilibration times were obtained from the 
non-steady-state version of this algorithm (refs 
B-16 and B-40; that was originally employed in the 
NASA studies The steady-state description per- 
mitted the use of a larger integration step size in the 
numerical algorithms and increased the solution 
speed significantly without decreasing the accuracy 
of the response for the long periods of time in 
which the investigators were interested 


Bom Marrow Controllar System 

An accurate description of the relationship 
governing erythropoietin release is not yet availa- 
ble, presumably because of the difficulty in measur- 
ing intrarenal oxygen levels and the uncertainty 
surrounding the specific location of the receptor 
cells The formulation used to relate tissue oxygen 
tension to erythropoietin release is in accord with 
the study of Adamson (ref. B-18), who found a 
semilogarithmic inverse relationship between daily 
urinary erythropoietin excretion and hematocrit in 
humans A parallel between urinary and plasma 
erythropoietin (ref B-7) as well as between 
hematocrit and tissue pOi (ref B-41) was assumed 
in deriving equation (B8) A similar relationship 
has been used by Hodgson (ref B-5). whereas Parer 


(ref B-42) has derived a linear relationship and 
Mylrea and AbbrecM (ref B-17) have used arterial 
oxygen-hemoglobin concentration (i t., Hb x S e O ) 
rather than tissue oxygen tension as the inde- 
pendent variable At the present time, sufficient 
data do not exist to reveal the precise shape of this 
function 

Measurement of plasma erythropoietin has. un- 
til recently, been restricted to levels above basal 
(ref B-l) Therefore, no data are available to con- 
firm the relationship to reduced release rates of 
erythropoietin This region is of particular interest 
because of its application of simulation to bed rest, 
space flight, and related disturbances in which 
chronic elevation of hematocrit follows plasma 
volume shifts. 

There is an abundance of information demon- 
strating that in experimental animals, a linear rela- 
tionship exists between red cell production and the 
log of erythropoietin concentration (eq (Bl 2) ) 
This has been observed, for example, in bioassay 
animals in which doses of erythropoietin are in- 
jected either singularly with iron uptake used as the 
index of ery thropoietic activity (refs. B-2 and B-24) 
or administered at frequent intervals for up to 
several weeks with production rate expressed in 
terms of increased red cell mass (refs B-9 and 
B-23) It is reasonable to assume that a similar dose- 
response relationship exists for the human, 
although confirmatory evidence is lacking In vivo 
estimates of the human function curve, especially 
for the suppressed erythropoiesis range, will be 
possible as erythropoietin becomes available in 
large quantities and as mote sensitive assay 
methods for this hormone are developed. 

The shape and position of the renal and bone 
marrow function curves have been found to be cru- 
cial elements in the control of erythropoiesi'. in 
general and in the long-term control of tissue ox- 
ygenation in particular Model parameters have 
been incorporated to allow for shifts in sensitivities 
and thresholds away from the normal operating 
points Values of controller sensitivities have not 
been well established in the human by direct 
methods and only to a limited extent in experimen- 
tal animils (ref B-5) Several studies suggest that 
alterations in these parameters occur during certain 
physiological stresses such as dehydration and hy- 
poxia (refs B-43 and B-44) and during pathological 
disturbances such as abnormal hemoglobin (ref 
B-34), eryihrocytosis (ref B-18). and hemolytic 
anemia (ref B-45) It appears that the rate of red 
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cell production is determined not onl> by the con- 
centration of erythropoietin but also by the size of 
the stem cell pool (ref B-46) If this is true, the 
bone marrow response to a given dose of 
erythropoietin should be greater than normal (i e.. 
an effective increase in Gi) The availability of iron 
to the erythron may also influence this function 
(ref. B-II) Certain hormones, such as androgens, 
as well as neural stimuli are assumed to exert their 
effect on erythropoiesis by their modification of 
erythropoietin release (i.e., an effective change in 
Gi.refs B- 12 and B-13). It is possible, using the cur- 
rent model, to predict these parameters within nar- 
row limits provided both the dynamic behavior of 
erythropoietin and red cell production ates are 
measured simultaneously during hematologic 
stress Unfortunately, such data are seldom availa- 
ble. especially for humans. 


Gantral Comments 

Feedback regulation of tissue oxy gen tension is 
accomplished solely by adjustments of hemoglobin 
lev els resulting from the output of a renal/bone- 
marrow controller Other parameters that are 
known tc effect acute changes in tissue oxygena- 
tion are incorporated explicitly in the model but are 
nonregulatory in nature and can be altered 
manually to test various hypotheses Similarly, the 
characteristics of the controller can also be adjusted 
to test their effect on long-term control of tissue 
pO: and red cell mass Such parameter adjustment 
(other than for the primary disturbance) has not 
been found to be essential, in most cases, to simu- 
late the basic behavior of the dynamic and steady- 
state response However, fine tuning of parameters 
is required to scale the model output and achieve 
closer agreement with experimental values In 
some cases, these studies indicate the need to pro- 
pose additional regulatory elements to provide, for 
example, a more realistic simulation of the 
erythropoietin response to hypoxia Other features 
of a general nature have also been identified that 
will increase the utility of the model even further, 
including (1) the effect of blood volume and 
viscosity on oxygen transport and (2) a description 
of stem cell kinetics and reticulocytosis (ref B-47) 


FORMULATION OF MOUSE MODEL 
System Parameter* 

The computer model for erythropoietic control 
was adapted to the mouse system by altering 
system parameters originally given tor the human 
to those that more realistically represent the mouse 
(ref B-48) Parameter values were obuined from a 
variety of literature sources as indicated in table 
B-II The immediate application of the mouse 
model was the study of the mouse as a potential ex- 
perimental model for space flight Data for the 
simulations were obtained from C. D R Dunn's 
experiments at the University of Tennessee 
Memorial Research Center and included studies o r 
dehydration and hypoxia The strain of mice used 
in these studies was C3H with an approximate 
weight of 25 grams Parameter values were chosen 
for this strain where possible. In certain cases, the 
literature values were superseded by values ob- 
tained directly from Dunn's studies In a few cases, 
mouse data were not available and data for the rat 
were substituted A comparison of system 
parameters for the mouse and human models is 
shown in tab’* B-II Aside from the obvious 
differences expected in fluid volumes, blood flows, 
and metabolic rates, larger differences were ob- 
served in the following erythrocyte lifespan (126 
days vs 20 days), 5 erythropoietin half-life (12 
hours vs 3.25 hours), and normal arterial pO: 
(12 67 kN/m : (95 mmHg) vs. 10 4 kN/m J (78 
mmHgi) The shorter lifespan of the mouse ted 
blood cells implies a turnover of erythrocytes 
which is sixfold faster That is, the daily rates of red 
cell production and destruction (as well as 
reticulocyte index) are approximately six times 
higher in the C3H mouse than the human * Other 
parameters which were found to be more similar 
between the two species were as follows hema- 
tocrit (40 vs 45). mean corpuscular hemoglobin 
concentration (0.375 vv 0.30). and maximum ox- 
ygen-carrying capacity of hemoglobin (1.34 vs 
1.41) 


’First and second numbers :n parentheses refer 10 human 
and mouse respectively 

•Typical values for mice red cel lifespan found in (he 
lueraiuri indicate only a threefold increase Values used here 
weie 'ounj b> Dunn 10 be much difTeren. in the C 3 H strain 
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TaBU B-ll — System Parameters for Erythropotesis Control Model 

Po'o^ftn btodc' Paromner oalui Rrlrrrm r Lniii 

im tr‘ 




Human 

Mo air 



Red cci: mass* 

RCM 

2000 

06 i 

B-49 

ml 

Plasma volume' 

n 

3000 

V 

B-49 

ml 

Blood volume 

BV 

5000 

1.40 

B-49 

ml 

Whole-body hcmatocni 

Ha 

400 

45 0 

B-49 

ml of picked RBC/100 ml 






t>kxx3 

Mean corpuscular hemoglobin 

HCHC 

.375 

300 

B-50 

g Hb/ml RBC 

concentration 1 






Hemoglobin concentration' 

Hb 

150 

12.5 

B- 50. B-51 

g Hb/100 ml blood 

0, capacity or blood* 

C.O 

201 

190 

B51 

ml CyiOO ml blood 

0, capacity of hemoglobin 

CHbO 

1 34 

1 41 

B-50. B-51 

ml 0./$ Hb 

pO, ten non at one-half Hb 

no 

27 

39 

B-52 

mmHg 

saturation' 


36 

52 


6 Pa 

•Arterial pO, 

P.0 

95 

71 

B 53 

minHg 



12 67 

104 


KPa 

Arterial Hb saturation 

5.0 

97 

99 

B-52 

Percent 

Renal metabolic rate* 


20 

04 

B 54 

ml O t 'min 

Renal blood flow* 

Q 

1200 

1 83 

B 55 

ml/min 

'Normal tissue pO- 

P.0 

20 

20 


mmHg 



2 67 

2 67 


kPa 

Eryihropoienn half-life* 

1 t-m 

12 

3 25 

B3b 

Hours 

Red cell lifespan' 

TRC 

126 

20 

B-17 

Days 

Erythrocyte maturation time* 

TBM 

4 

35 


Days 

Normal RBC production rate* 

RCP 

22 

0437 

B-49 

ml RBC/day 

RBC turnover rate 

RKC 

1 1 

6 93 

B-49 

Pcrccnt/day 
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Although ihc *rteri«! pO: in the moose is much Vglues for renal blood flow of the mouse were 

lower than in the human, the oxygen saturation of not available, and data from rats were used (6 ml 

hemoglobin of both species is nearly identical (97 rntn“' g~' tissue) 

percent vs 99 percent) This is a result of the dis- 
tinctly different oxygen-hemoglobin dissociation 

curves shown in figure B-8 and reflecied in the .Scaled Parameters 

different PS 0 values (3 53 kN/m : (26.5 mmHg) vs 

5.2 kN/m 1 (39.0 mmHg)) The P50 differences im- Some parameters of the mouse model differ con- 

ply that, at the same level of tissue oxygen tension, siderably from the human model because of scaling 

oxygen is more easily unloaded in the mouse than factors alone The values used in the model are 

in the human I: should be noted that the normal given on an absolute basis for the whole animal 

pOi of arterial blood assumed here (10 4 kN/m ; (78 rather than as a specific propert} in terms of ' per 

mmHg)) was obtained from rat data (refs B-51 and gram of tissue" In terms of specific units, the 

B-56) and has been used in a previous model vali- differences between the mouse and human system 

dated for the mouse with reasonabl) good results are much smaller, as shown m table B-lll 

(ref. B- 1 7) No corresponding mouse data could be 
located 
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FIGURE R-6 — 0*r|rn-b«B>o*lobin equilibrium runes f or human ut mouse tapeilmental versus theoretical values 


Table B-lll.— Absolute Venus Spec if n Parameter Values 


Parameter 


Absolute units 


Sptetfk units * 


Human 

Mouse 

Units 

Human 

Mouse 

Units 

Rer 1 oel' mast 

2000 

06) 

ml 

lib 

252 

mlA| ^od> *t 

Ft* » na volume 

MOO 

77 

ml 

42 9 

Ml 

mlA| bod) »t 

P jod volume 

S00U 

1 40 

ml 

71 4 

560 

mlAg bod) *t 

llenai blood flow 

1200 

1 1) 

ml/mm 

4 21 

6 10 

ml mm* 1 1 * 1 

tissue 

Renal 0, consumption 

20 

04 

mt/min 

07) 

D) 

ml mm* ' g~ 1 

i issue 

Rod> 0, oonaumpnon 

250 

51 

ml/mm 

00)5“* 

0255 

m! mm" ' g" 1 

tissue 
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Oxygen Balance 

The balance of oxygen supply versus oxygen de- 
mand is crucial to the feedback regulation of 
erythropoiesis A parameter reflecting this com- 
plex balance is the tissue oxygen tension that is 
believed to govern the release of erythropoietin 
The oxygen balances for the human and mouse 
systems as used in the model are given in table 
B-IV. 

Oxygen consumption per gram of renai tissue in 
the mouse is approximately twice that for the 
human (Overall total oxygen consumption per 
gram body weight is nearly seven times greater in 
the mouse.) This higher oxygen demand of the 
mouse is satisfied in two ways in the model First, 
there is a 50 percent greater efficiency in oxygen 
extraction as indicated in table B-IV. (Note that, in 
both species, the amount of oxygen delivered at 
rest is more than sufficient; i.e . roughly 10 times 
that required by the tissues ) Second, there is a 30 
percent higher blood oxygen supply per gram of 
tissue because of greater tissue blood flow in the 
mouse 


Table B-l\ — Oxygen Balance ai Kidney 


0 . supp'y equals 
BF * 0 :oni enaction 


Kidnf , 


0- demand 


(X venous 


Parameter 

Human 

Mouse 

L'nns 

Oiygen demand 

70 

004 

ml Oj-'mtn 


O’J 

IS3 

ml 0, min' ' g' 1 

Otyger. supply 

pO. arterial 

95 

71 

mmHg 


12 6’ 

10 4 

kh 

sO. inernl 

97 4 

9|fr 

Percent saturation 

0 : eoncentraiior 

I9t 

III 

ml O/liter blood 

BF 

1200 

1 13 

ml Mood'min 

0, supply rait 

235 

34 j 

ml 0,/mm 


139 

1 143 

ml P,min''g"' 

Osygen venous 

pO. venous 

56 

57 

mmHg 


7 4' 

’6 

kPa 

lO,. venous 

K 

<6 

Percent saturation 

pO. tissue 

20 

20 

mmHg 


2 6' 

7 6' 

kPa 

Percent o»>ger 

1 ' 

134 

Percent 

eitramor • t(i 
derrundt 'n' surr'< 1 


The normal tissue oxygen tension is arbitrarily 
assumed to be identical in both model systems, i.e , 
2 67 kN/m ; (20 mnHg) The equation describing 
oxygen diffusivity to the tissues from venous 
capillaries is given in the steady state as 

Net oxygen delivery * tissue oxygen consumption 

= ^p0 2 ,vein - pCL, tissue^ X K 

where K - conductivity coefficient - Oj diffu- 
sivity times the capillary surface area The ratio 
^min^mouse WOu, d be expected to reflect the sur- 
face area ratio between species if diffusivity is 
assumed c : milar in mouse and man Therefore, if 5 
is capillary surface area, then 


nioufc mou* 

_ (O; wUKv-imi.hon)^ ^ (PO; ™ pO ; 

‘ (O- con$u.npi,.in) mouw (p0 2 vein pO. tissue)^, 


70 ml mm (57 5 20) nimHt 

OCMmi'mm (St* 20lmmH^ 

* MS 


This is in good agreement with the surface area 
ratio of 650 of the glomerulus derived from data in 
reference B-54 (p P4), lending support to the 
general representation of the kidney in the com- 
puter model. 


Functional Relationships 

Three functional relationships are included in 
the computer model (1) oxygen-hemoglobin 
equilibrium curve (OEC), (2) erythropoietin 
release as a function of tissue pOi, and (3) 
erythrocyte production rate as a function of 
erythropoietin concentration The first of these is 
shown in figure B-8 and will subsequently be de- 
scribed in detail The form of the function curves 
for erythropoietin and red cell release will be 
assumed identical in the mouse and human 
models There is no reason at the present time to 
take issue with this assumption, particularly since 
the bone marrow function was originally obtained 
from the mouse These curves (as shown in figs 


B-15 




C'.g- if' rVisi- < 
OF tOCH QUALITY 


B-2 and B-3 and as used in the models) are repre- 
senied in normalized form (i c . percent of control) 
so that any species may be represented The gain 
factors G and G. representing the slope of the rela- 
tionships, may be different between species This is 
of little concern in the basic design of the model 
because these parameters will be adjusted during 
the simulation process and their actual values will 
be estimated by ‘Titling" the model output to the 
experimental data 

The equation describing the sigmoidal OEC is a 
form of the Hill equation and is shown in the insert 
of figure B-8 The two solid lines represent human 
and mouse blood, respectively, and were recently 
obtained from blood samples of the C3H mouse 
(ref B-49) The value of /*50 is explicitly stated in 
the equation so that shifts in oxygen-hemoglobin 
affinity may be easily described The value of the 
exponent k, found from the best fit of the mouse 
curve, also provides a good fit of the human curve 
as shown in figure B-8 Thus, the only difference 
between the equation describing the human and 
mouse OEC is the value of PSO. 
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